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Abstract. We present new broad-band optical and near-infrared CCD imaging together with deep optical 
intermediate-resolution spectroscopy of Mkn 1087 and its surrounding objects. We analyze the morphology and 
colors of the stellar populations of the brightest objects, some of them star-formation areas, as well as the kine- 
matics, physical conditions and chemical composition of the ionized gas associated with them. Mkn 1087 does not 
host an Active Galactic Nucleus, but it could be a Luminous Compact Blue Galaxy. Although it was classified as 
a suspected Wolf-Rayet galaxy, we do not detect the spectral features of these sort of massive stars. Mkn 1087 
shows morphological and kinematical features that can be explained assuming that it is in interaction with two 
nearby galaxies: the bright KPG 103a and a dwarf {Mb ~ —18) star-forming companion. We argue that this dwarf 
companion is not a tidal object but an external galaxy because of its low metallicity [12-|-log(0/H) — 8.24] with 
respect to the one derived for Mkn 1087 [12-|-log(0/H) = 8.57] and its kinematics. Some of the non-stellar objects 
surrounding Mkn 1087 are connected by bridges of matter with the main body, host star-formation events and 
show similar abundances despite their different angular distances. These facts, together their kinematics, suggest 
that they are tidal dwarf galaxies formed from material stripped from Mkn 1087. A bright star-forming region at 
the south of Mkn 1087 (knot #7) does not show indications of being a tidal galaxy or the product of a merging 
process as suggested in previous works. We argue that Mkn 1087 and its surroundings should be considered a 
group of galaxies. 



Key words, galaxies: starburst — galaxies: interactions — 
galaxies: abundances — galaxies: kinematics and dynam- 
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1. Introduction 

One of the most controversial issues about Hii galax- 
ies is the triggering mechanism of the violent bursts of 
star formation that are observed in them. The estimated 
star formation rates in these starburst galaxies are so 
high that the material available for the creation of new 
stars would be exhausted very soon compared to the 
age of the uni verse. T his problem is more marked in 
dwarf galaxies. iThuan I (jT991) suggested that the bursts 
have an intermittent behavior. However, star formation 
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Newton Group of Telescopes and Instituto de Astrofi'sica 
de Canarias in the Spanish Observatories of Roque de Los 
Muchachos and Teide of the Instituto de Astrofi'sica de 
Canarias. 



episodes show a l arge variability in duratio n, ranging from 
10^ - 10^ vear llRieke fc Lebofskl Il985^ to more than 
10^ yr l|Hunter fc GallagheJll985|) . In fact, the age of the 
bursts and the galaxy that hosts them is also a controver- 
sial problem. The knowledge of both the triggering mech- 
anism and the age of the starbursts would help to under- 
stand the evolution of the galaxies. 



Hierarchical forma tion models of galaxies 
llKauffmann et all Il997l) predict that most galaxies 
have formed by merging of small clouds of protogalactic 
gas. In this sense, studies of metal deficient blue compact 
dwarf galaxies (BCDGs) embedded in large presumably 
primordial Hi clouds provide excell ent clues about the 
trigger and youth of the galaxies l|Thuan et al.l Il995t 
llzotov fc ThuanI Il999l) although most of them show 
an u nderlying old stellar component ijNoeske et alJ 
l200.'^ . Some authors ijSanders et al.lll98^ proposed that 
galaxy interactions are the main triggering mechanism 
in the starbursts of dwarf and irregular galaxies. The 
interactions are not usual with nearby giant galaxie s 
( Camoos-Aguilar et alJ Il99.'4 iTelles fc TerlevichI Il995l) . 
but they seem to be more frequent with low surface 
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Table 1. Summary of observations 



Observations 


Telescope 


Date 


Exp. Time 


Spatial 


Filter/ 


P.A. 


Spectral 


AA 








(s) 


( Pix ) 


grating 


( ) 


(A pix ) 


(A) 


Broad-band 


2.2m CAHA 


00/12/19 


3 X 1200 


0.53 


V 








imaging 


2.56m NOT 


03/01/20 


3 X 300 


0.19 


B 










2.56m NOT 


03/01/20 


6 X 300 


0.19 


R 










1.55m GST 


02/09/24 


120 X 20 


1.0 


J 










1.55m GST 


02/09/24 


240 X 10 


1.0 


H 










1.55m Gb 1 


02/09/24 


360 X 5 


1.0 


Ks 








Intermediate 


4.2m WHT 


00/12/31 


1800 X 3 


0.20 


R600B 


358.0 


0.45 


3650-5100 


resolution 


4.2m WHT 


00/12/31 


1800 X 3 


0.36 


R136R 


358.0 


1.49 


5300-6650 


spectroscopy 


2.5m INT 


99/12/28 


1200 X 4 


0.40 


R400V 


357.0 


1.40 


3500-7700 




2.5m INT 


99/12/29 


1200 X 3 


0.40 


R400V 


6.0 


1.40 


3500-7700 




2.5m INT 


99/12/29 


1200 X 3 


0.40 


R400V 


62.0 


1.40 


3500-7700 




2.5m INT 


99/12/29 


1200 X 3 


0.40 


R400V 


120.0 


1.40 


3500-7700 



brigh tness galaxies I Tavlor et al.1 Il99fit IWilcots et alJ 
Il99fil) . 

On the other hand, the gas-rich dwarf irregular galax- 
ies could have formed by material ejected from the disks 
of the parent galaxies to the interga lactic medium by tidal 
forces l|Qkazaki fc Taniguchillioootl . The recycled objects 
formed in this way (the tidal dwarf galaxies, TDGs) give 
import ant clues to both ag e and past evolution of the 
galaxy JDuc fc Mirabe]lll998l) . In some cases, star-forming 
TDGs are foun d at the end of long st ellar tails in inter- 
acting systems llWeilbacher et al.|l2003|) . 

In one way or in another, the interactions seem to play 
a key role in the triggering of the star formation in dwarf 
galaxies. An increasing number of starburst galaxies show 
Wolf-Rayet (WR) features in their integrated spectra. The 
so-called Wolf-Rayet galaxies have a broad emission fea- 
ture around 4650 A that usually consist of a blend of lines 
He II A4686 (the most important line), Ciii/Civ A4650 
and Niii A4640, and it indicates the presence of a sub- 
stantial population of this sort of massive stars. The most 
massive OB stars became WR stars around 2 and 3 Myr 
since t heir birth, and spend bet ween 3 and 6 Myr in this 
phase l|Maeder fc Mevnetlll99^ . The short time in which 
the WR features can be observed makes these systems 
ideal objects to look for the probable interaction mecha- 
nism of the star formation. They also provide information 
about the age and prop erties of the starburs ts. Studying a 
sample of WR galaxies. iMendez fc Esteban (,200QJ (here- 
inafter MEOO) suggested that interactions with or between 
dwarf objects could be the main star formation triggering 
mechanism in dwarf galaxies. Once triggered, the evolu- 
tion of a starburst depends on the details of its interaction 
with the galaxy environment. MEOO also remarked that 
the interacting and/or merging nature of the WR galax- 
ies can only be detected when deep and high-resolution 
images and spectra are available. 

Mkn 1087 (II Zw 23, UGC 3179) was firstly described 
by IZwickvl l)l97l|) as an elliptically shaped BCG with a 
very extended filament to the north and a less exte nded 
one to the south. These filaments were studied by iKeell 



119881) using optical imaging and spectroscopy and radio 
continuum mapping of the system. He found that the fil- 
aments were redshifted with respect to the central zones 
of the galaxy, suggesting that they are the result of neu- 
tral gas infalling from a very massive Hi cloud in which 
Mkn 1087 is embedded. He proposed that Mkn 1087 is a 
galaxy in formation and suggested that the star formation 
process has been occurring continuously for the past 10^ 
year. He also remarked that it could be as a consequence 
of the interaction with a companion galaxy located 10" 
to the south of the Mkn 1087 nucleus (our knot #7, see 
§3.1). Mkn 1087, which is at around 111 Mpc (assuming 
Hq= 75 km s~^ Mpc at this distance 1"=538 pc), has 
a companion galaxy named KPG 1 03a, at 81 kpc to th e 
southwest and at th e same redshift ijMarzke et al.lll99^ . 
iMalkan et all l|l998() obtained a HST WFPC2 image of 
Mkn 1087 showing that it is a spiral galaxy with dust 
lanes on one side of its major axis. 

The presence of WR features in the integrated 
spectrum of Mkn 10 87 was first investigated by 
iKunth fc Jouberd l|l98 5V who gave an upper limit to 
the flux of th e blue WR bump around Hen A4686. 
IVacca fc Contil l|l992|) did not measure the Heii A4686 
emission line but reported an upper limit of 1600 to the 
total late- type nitrogen-rich WR stars (WNL ) cont ained 
in the central region of Mkn 1087. IVaceh et all l|l997l) per- 
formed a spectroscopical study of 46 galaxies with emis- 
sion line spectra in where Mkn 1087 was included, and 
reported a detection of the He ii A4686 emission line, giv- 
ing its flux with an error less than 50%. However, they 
did not give information about its line width, so the neb- 
ular or stellar origin c an not be specified. As consequence, 
ISch aerer et alJ l)l999|) classified Mkn 1087 as a suspected 
Wolf-Rayet galaxy. 

MEOO presented deep CCD UBV images of Mkn 1087, 
finding and characterizing several faint objects connected 
by bridges and tails with the main bo dy of the galaxy. 
They used ILeitherer fc HeckmanI l|l995|) population syn- 
thesis models to estimate the ages of the emission zones 
detected in their images, finding that some of these emis- 
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Fig. 1. Deep optical image of Mkn 1087 and its surrounding in R filter. Mkn 1087, KPG 103a, the north companion galaxy 
(N) and the rest of the surrounding objects are labeled. The background cluster of galaxies is indicated with a circle. 



sion knots are around 5 Myr old. MEOO interpreted Mkn 
1087 as a galaxy experiencing an interaction event with 
one or more of the neighboring objects. The main aim of 
this study is to perform a more complete observational 
analysis of these surrounding objects in order to establish 
their nature. 

2. Observations and data reduction 

2.1. Optical imaging 

Optical V images of Mkn 1087 were carried out on 2000 
December 19 with CAFOS (Calar Alto Faint Object 
Spectrograph) in image mode at the Cassegrain focus of 
the 2.2m telescope of the Centre Astronomico Hispano- 
Aleman (CAHA) at Calar Alto Observatory (Almeria, 
Spain). A SITe CCD detector (2048 x 2048 pixels) with a 
pixel size of 24 /im and spatial resolution of 0.53" pixel^^ 
was used, as well as the standard Johnson V filter. Three 
1200 seconds exposures were added for each filter to ob- 
tain a good signal-to-noise and an appropriate removal of 
cosmic rays in the final images. The measured full width 
at half maximum (FWHM) of the point-spread function 
(PSF) was approximately 1". Images were taken under 



photometric co nditions, and the standard field 98-1119 of 
iLandoltl (^Q^) was used to perform the flux calibration. 

Images in B and R filters were taken on 2004 January 
20 at the 2.56m Nordic Optical Telescope (NOT) at 
Roque de los Muchachos Observatory (La Palma, Canary 
Islands, Spain). We used the ALFOSC (Andalucia Faint 
Object Spectrograph and Camera) in image mode with a 
Loral/Lesser CCD detector (2048 x 2048 pixels) with a 
pixel size of 15 /im and spatial resolution of 0.19" pixel"^, 
and the standard Johnson filters. We obtained three 300 
seconds exposures for B filter and six 300 seconds expo- 
sures for R filter, that were combined to obtain the final 
images with F WHM of the P SF of 0.8". The standard 
star SA97-284 l|Landolijll99l|) was used to flux calibrate 
the final images. 

Twilight images of different zones of the sky (blank 
fields) were taken for each filter in order to perform the 
flat-field correction. All the reduction process (bias sub- 
traction, fiat-fielding and flux calibration) was done with 
IRAF^ package following standard procedures. IRAF soft- 
ware was also used both to determine the 3(t contours 



^ IRAF is distributed by NOAO which is operated by AURA 
Inc., under cooperative agreement with NSF 
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over the sky background level and to obtain the photo- 
metric values of the galaxies and the knots surrounding 
Mkn 1087. 

We have also re-analysed the previous U, B and V 
images obtained by MEOO at 2.56m NOT (Nordic Optical 
Telescope) on 1997 February 6. See MEOO for a detailed 
description of these observations. 

2.2. Near-Infrared imaging 

We used the 1.55m Carlos Sanchez Telescope (CST) at 
Teide Observatory (Tenerife, Canary Islands, Spain) to 
obtain the near-infrared images on 2002 September 24. 
We used CAIN camera (256 x 256 pixels) with a pixel 
size of 40 /im and a spatial resolution of 1" pixel^^ in the 
wide mode to obtain images in J (1.2 /im), H (1.6 /im) 
and Ks (2.18 /im) broad-band filters. We took 10 series of 
6 consecutive individual 20 s exposures in J filter, 10 series 
of 12 individual 10 s exposures in H filter and 10 series of 
12 exposures of 5 s duration in Ks . Each series of exposures 
was made at slightly different positions to obtain a clean 
sky image. In this way, we obtained a single image of 20 
minutes in J and H filters, and a single one of 10 minutes 
in Kg filter. We repeated this procedure two times for J 
and H and three times for Kg, combining these images in a 
definitive one for each filter to increase the signal-to-noise 
ratio. The FWHM of the PSFs of these final images was 
2.8", 2.6" and 2.5" for J, H and K^, respectively. Bright 
and dark dome flat-field exposures were taken for each 
filter that were combined to obtain a good flat-field im age. 
The standard stars AsOl and Asll l|Hunt et al.ll998^ were 
used for flux calibrate the images. 



2.3. Intermediate-resolution spectroscopy 

We have obtained long-slit intermediate-resolution spec- 
troscopy in five slit positions centered on the nucleus of the 
galaxy and covering most of the surrounding objects. We 
have used the 4.2m William Herschel Telescope (WHT) 
and the 2.5m Isaac Newton Telescope (INT), both at 
Roque de los Muchachos Observatory (La Palma, Canary 
Islands). 

INT spectroscopy was carried out on 2000 December 
28 and 29 with the IDS spectrograph at the Cassegrain 
focus with the 235 mm camera. An EEV CCD with a 
configuration of 2148 x 4200 pixels of 13.5 /im was used. 
The slit was 2.8' long and 1" wide. The grating R400V 
was used, with a dispersion of 104.5 A mm~^ and effec- 
tive spectral resolution of 1.40 A. Each individual spec- 
trum covers from 3200 to 7700 A. The spatial resolution 
is 0.40" pixel"^. Four slit positions of Mkn 1087 were ob- 
served at different position angles, which were chosen in 
order to cover different morphological zones. Three or four 
20 minutes exposures were taken and combined to obtain 
good signal-to-noise and an appropriate removal of cosmic 
rays in the final spectra for each slit position. Comparison 
lamp exposures of CuAr were taken after each set of spec- 



tra. The absolute flux calibration was achieved by obser- 
vations of the standard stars Feige 56, Hiltner 600 and 
Feige 110 llMa,ssev et a,].lll98Sh . 

An additional slit position of Mkn 1087 was obtained 
using 4.2m WHT. It was carried out on 2000 December 
31 with the ISIS spectrograph at the Cassegrain focus. 
Two different CCDs were used at the blue and red arms 
of the spectrograph: an EEV CCD with a configuration of 
4096 X 2048 pixels of 13 /xm in the blue arm and a TEK 
with 1024 X 1024 of 24 /im in the red arm. The dichroic 
used to separate the blue and red beams was set at 5400 
A. The slit was 3.7' long and 1" wide. Two gratings were 
used, the R600B in the blue arm and the R316R in the red 
arm. These gratings give reciprocal dispersion of 33 and 
66 A mm , and effective spectral resolutions of 2.0 and 
3.9 A for the blue and red arms, respectively. The blue 
spectra cover from 3600 to 5200 A and the red ones from 
5400 to 6800 A. The spatial resolution is 0.20" pixer^ in 
the blue and 0.35" pixel"^ in the red. We obtained three 
30 minutes exposures in each arm that were combined to 
obtain the final blue and red spectra for this slit posi- 
tion. Comparison lamp exposures of CuAr for the blue 
arm and CuNe for the red one were taken after each set 
of spectra. The absolute flux calibration was achieved by 
obs ervations of the sta ndard stars G191, B2B and Feige 
34 (jMassev et al.lll98'8l) . WHT observations were made at 
air masses around 1.1. INT observations were made at air 
masses between 1.1 and 1.3. Consequently, no correction 
was made for atmospheric differential refraction. 

All the CCD frames for spectroscopy were reduced us- 
ing standard IRAF procedures to perform bias correc- 
tion, flat-fielding, cosmic-ray rejection, wavelength and 
flux calibration, and sky subtraction. The correction for 
atmospheric extinction was performed using an average 
curve for the continuous atmospheric extinction at Roque 
de los Muchachos Observatory. For each two-dimensional 
spectra several apertures were defined along the spatial 
direction to extract the final one-dimensional spectra of 
each knot, centering each aperture at its brightest point. 
For WHT observations, identical apertures were defined 
on the "blue" and "red" frames. Small two-dimensional 
distortions were corrected fitting the maxima of [O ii] 
AA3726,3729 doublet and Ha in INT frames, using also 
H/3 and continuum emission for WHT ones. 

The journal of all the imaging and spectroscopical ob- 
servations can be found in Tabled 



3. Results 

3.1. Optical and near-infrared imaging 

T he deep CAH A V image of the system was presented 
in lEsteban et a l. (200^, showing the main non-stellar ob- 
jects surrounding Mkn 1087 found by MEOO. We desig- 
nated the objects following those authors. Several bridges 
can be observed connecting the main body of Mkn 1087 
with the objects #1, #2, #3, #11 and #12, the last two 
ones located in the direction of the nearby galaxy KPG 
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Table 2. Results of optical and near-infrared aperture photometry in Mkn 1087 and KPG 103a knots. 
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E{B - V) 








{U-B)o 


{B - V)o 


(V - R)o 


(V - J)o 


{J-//)o 


{H-Ks)o 


iViK.Il J.UO 1 


ly^ 


13 78 




1 9 qQ_|_n fiA 


n zLi -l-fi ns'' 




n 9n-l-n ns 


n c;9-i-n rifi 

U. OZ ^U. UD 


n 9n-i-n riR 




c 


0. 14 


15.19 


i4.DZ±U.Uo 


14. Jo±U.Uo 


— U.57±u.uD 


— O.U2±O.Uo 


nil _l_ r\ r\£2 
U. ii±U.UD 


U.oi±U. iU 


U.io±U. iU 


n 1 c _i_ n 1 n 

U.i5±U.iu 


N 


0.10 


17.58 




1 T n 1 _Ln nt; 
i / .Uiihu.UD 




— U.UD±U.UD 


U.i4±U.iU 


n n 1 _Ln no 


n 1 _Ln no 
U.i8±U.U8 


n 1 _Ln no 
U. io±U.U8 


#1 al 


0.07^= 


19.18 


18.89±0.05 


>18.7 


— U.7d±u. Id 


— ().()1±U.U8 


U.10±0.U8 


<0.2 






#2 


0.07"= 


20.01 


19.72±0.05 


>18.7 


-0.78±0.15'' 


0.05±0.08 


0.12±0.08 


<1 






#3 b2 


0.07= 


18.32 


18.03±0.03 


17.28±0.06 


0.08±0.30'' 


0.11±0.06 


0.26±0.06 


0.64±0.10 


0.5±0.2 




#7 W 


0.16 


18.29 


17.63±0.04 




-0.87±0.08'' 


-0.04±0.08 


0.03±0.08 








#8 a4 


0.08 


19.10 


18.77±0.05 




— 0.27±0. 12'' 


0.07±0.10 


0.17±0.08 








#11 


0.07= 


20.69 


20.40±0.05 






0.35±0.10 


0.08±0.10 








#12 c3 


0.07= 


19.61 


19.32±0.05 


18.5±0.1 


-0.1±0.3 


0.31±0.09 


0.24±0.08 


0.55±0.15 


0.4±0.2 




#13 a2 


0.40 


19.25 


17.61±0.05 




-0.64±0.12 


0.01±0.09 


0.04±0.08 








#14 cl 


0.07 


18.45 


18.16±0.05 




-0.59±0.12 


0.28±0.10 


0.23±0.10 








#15 d2 


0.19 


20.90 


20.12±0.08 




-0.6±0.3 


0.22±0.12 


0.23±0.10 








KPG 103a 

Kl 

K2 

K3 


0.07= 
0.07= 
0.07= 
0.07= 


14.05 
17.55 
17.10 
17.41 


13.77±0.04 
17.26±0.05 
16.81±0.04 
17.12±0.04 




0.17±0.08'' 
-0.40±0.10 
-0.49±0.08 
-0.39±0.08 


0.18±0.06 
0.18±0.10 
0.12±0.06 
0.19±0.06 


0.27±0.06 
0.20±0.08 
0.11±0.06 
0.33±0.06 









° Average of the values obtained for C, #7, #8, #13, #14 and #15. 
^ Adopting the mu values from M EOO. 

= Galaetic reddening adopted from lSchlegel et al.l 1199811 . E(B - V) = 0.07. 




103a. Knot #7 seems to be an intense star-formation zone 
off-center Mkn 1087. 

In order to obtain a deeper image, we carried out new 
observations of the system in R broad-band filter. We 
present this image in Figure ^ The bridges commented 
before are now clearly detected, especially the ones that 
end at #11 and #12. The bright bridge connecting Mkn 
1087 with #1 is apparently splitted in several small knots. 
The north companion object seems to be clearly elliptical 
and shows a very weak plume in the direction to knot #1. 
Although bridges between the two brightest galaxies are 
not detected, a long very weak plume can be observed at 
the east of KPG 103a. Another diffuse long plume is de- 
tected at the west. It is important to note the high number 
of (apparently) background non-stellar objects (presum- 
ably galaxies) that are detected in the new deep R image 
(we have estimated around 60 candidates in the entire fi- 
nal R frame). They are uniformly distributed in the field 
of view of the image (about 6') and are most probably 
background galaxies. We also want to remark the cluster 
of galaxies found at the upper right corner of the image, 
with central position around a = 04h 49min 35sec, S = 
+03° 21' 04". It is not catalogued in the NASA/IPAC 
Extragalactic Database (NED). We analyze it in §4.7. 

Some important structures inside the main galaxies are 
now clearly distinguished. In Figure [3 we show the cen- 
tral part of Mkn 1087 (up) and KPG 103a (down). The 
center of Mkn 1087, C (knot #5 following MEOO), is the 
brightest knot of the system. Mkn 1087 reveals two spi- 
ral arms or structures: an upper arm connecting with the 
bridge towards knot #2 and a bottom arm that seems to 
end in knot #7. Furthermore, some new knots are now de- 
tected inside the main body of both galaxies. Especially 
intense arc the three bright knots Kl, K2 and K3, that are 
resolved in the center of KPG 103a. They show blue col- 
ors (see Table EJ that are representative of young popula- 



Fig. 2. Detailed images in R band of Mkn 1087 (up) and KPG 
103a (bottom). The brightest knots inside both galaxies are 
labeled, as well as the external objects #2, #3 and #12. A 
weak plume at the east of KPG 103a, that curves toward the 
south, is also barely observed (it is indicated by an arrow). 
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A R.A. 

Fig. 3. Logarithmic contour maps in B and J filters of Mkn 1087. M indicates the main body of Mkn 1087, C the center of it, 
K the nearby galaxy KPG 103a and N the north companion galaxy. The contours are 3a, 10a, 20a, 50a and lOOa over the sky 
level for both maps. The slit positions used for spectroscopy are also shown over the B image, indicating with ticks the center 
of the different zones from which we have extracted the one-dimensional spectra. 



tion, suggesting that these knots are star-formation areas 
in the nucleus of KPG 103a. The WFPC2 H ST image of 
Mkn 1087 presented bv'Malka n et alJ l)l998l) clearly shows 
knots C, #7, #8, #13 and #14 and the arms of the spiral 
disk we have distinguished in our NOT images. The HST 
image suggests that object #2 could be a background ob- 
ject that is coincident with the eastern bridge of Mkn 1087, 
although this cannot be confirmed without a spectrum 
of the object. However, this knot shows very blue colors, 
(U-B) = -0.78, similar to other knots around Mkn 1087, 
so we prefer to consider knot #2 in our study. 

In Figure 01 we show the logarithmic contour map of 
Mkn 1087 and its surroundings in filters B and J. We 
clearly detect the north companion object in all NIR im- 
ages. However, the knots #3 and #12 are only detected 
in J and H filters. Knots #1 and #2 are barely detected 
in the J image. We did not observe the galaxy KPG 103a 
in NIR. 

MEOO performed aperture photometry of the knots of 
the system in U, B and V filters. We have improved their 
analysis performing the aperture photometry of the com- 
panion object in our new B, V and R images. We have 
determined the area for which we have integrated the flux 
using the 3a level isophote over the average sky level in 
the B image. The same integration area was then used 
for V and B filters. We have also measured the rest of 
the knots in the system, finding similar values to those 
obtained by MEOO. In addition, we have determined the 
photometric U, B and V values of some new interesting 



areas (#13, #14 and #15) because they correspond to lo- 
calized bright nebula emission in the disk of Mkn 1087 as 
our spectroscopy reveals (see next section) , they probably 
correspond to giant Hii regions. Finally, we have analyzed 
the NIR images to obtain the photometric J, H and A'^ 
values of the detected objects. The integration area in J 
filter (defined as the ia level isophote over the sky level) 
was used for each NIR filter. We also estimated a lower 
limit for to,/ in knots #1 and #2. All these new opti- 
cal and NIR data are shown in Table |21 We only give 
upper limits to the J mag nitudes for knots #1 and #2. 
iHunter fc Gallagheil lIlDSfij) performed JHK photometry 
within 23" aperture and found J — H = 0.65, H — K ^ 
0.28 and V — J = 1.76, suggesting a dominant contribution 
of old populations at these wavelengths. 

The reddening constant, C(H/3) — obtained from the 
absorption-corrected intensities of H i Balmer lines in our 
optical spectra of each member — was used to correct 
the photometric data for the extinction by interstellar 
dust. We describe the method used to obtain C(H/3) from 
the spectra in the ne xt section. The relation obtained by 
iKaler fc Lut j l| 19851) between C(H/3) and the extinction 
in V filter {Ay) was used [Ay ^ 2C(H/3)], and a stan- 
dard ratio of Ry = Ay I E{B — V) — 3.1 was assumed. 
Calzetti et at~l l|200 0l) propose a starburst reddening curve 
with Ry = 4.05 ±0.80. If we use this curve for our objects, 
the B magnitudes would be between 0.08 and 0.24 brighter 
than the ones we have calculated assuming Ry = 3.1. 
We would also obtain bluer colors, but the differences are 
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Knot 


C WHT 


C INT 


N WHT 


N INT 


7 WHT 


#7 INT 


#1 


#3 


#8 


#13 


#14 


#15 


-Mb 


20.61 


20.61 


18.06 


18.06 


17.60 


17.60 


16.34 


17.20 


16.46 


17.62 


17.07 


15.11 


± 


0.03 


0.03 


0.03 


0.03 


0.04 


0.04 


0.05 


0.03 


0.05 


0.05 


0.05 


0.08 


Slit ( X 1 ) 


5.4 


12 


3.6 


8.4 


6.0 


3.6 


5.6 


4.0 


5.2 


7.2 


5.6 


8.0 


Distance (")^ 








77.0 


77.0 


9.8 


10.0 


47.6 


35.2 


10.9 


10.8 


8.0 


19.2 




0.17 


0.25 


0.17 


0.12 


0.24 


0.25 


0.10 


0. 10 


0.12 


0.62 


0.11 


0.29 


± 


0.02 


0.03 


0.03 


0.03 


0.03 


0.03 


0.03 


0.03 


0.03 


0.12 


0.04 


0.05 




1.7 


2.4 


0.2 


0.9 


2.9 


1.5 






1.1 


2 


5.6 


2 


± 


0.1 


0.3 


0.2 


0.2 


0.5 


0.3 






0.3 




1.7 




Te [O III]*^ 


0.71 


0.74 


1.09 


1 13 


77 


78 


82 


76 


88 


82 


86 


81 


tI[o II]° 


0.80 


0.82 


1.06 


1.09 


0.84 


0.84 


0.88 


0.83 


0.92 


0.87 


0.90 


0.87 




220 


170 


115 


150 


120 


210 






150 


120 


<100 


<100 










+ 117 


+ 110 


+ 118 


+ 122 


-63 


+8 


+ 125 


-85 


-70 


-75 


± 






15 


27 


15 


27 


27 


80 


27 


27 


35 


35 


F(H/3)'5 


133 


184 


4.58 


11.3 


16.6 


1.76 


3.05: 


1.58: 


6.03 


1.92 


8.46 


0.94: 


± 


4 


6 


0.16 


0.6 


0.7 


0.09 






0.37 


0.26 


0.62 






87.4 


87.3 


53.6 


66.5 


330 


256 


110 


20: 


72.0 


30.9 


140 


39.7 


± 


2.6 


3.4 


2.1 


8.2 


25 


28 


25 




4.5 


2.8 


20 


9.5 


W(H/3)° 


22.3 


22.9 


25.0 


27.4 


62.3 


66 






16.6 


5.6 


52 


8.0: 


± 


0.9 


2.1 


1.6 


2.4 


5.4 


12 






2.6 


1.3 


8 




W(H7)'^ 


8.1 


6.6 


8.2 


5.3 


34.3 


20.2 






8.3 




18 




± 


0.4 


0.6 


0.9 


0.5 


4.9 


2.2 






1.6 




5 




w{[o m]r 


16.7 


20.5 


144 


46.6 


116 


63.2 






25.7 


6.8 


42.0 


7.8: 


± 


0.6 


0.9 


18 


3.7 


9 


3.6 






2.8 


1.2 


3.6 




Ago" 


6.0 


6.0 


7 


6 


5.0 


4.5 


6 


150 


6 


8.5 


5 


8.0 


± 


0.5 


0.5 


1 


1 


0.5 


0.5 


2 


30 


1 


1.5 


1 


1.5 



^ Distance from the center of Mkn 1087. 

^ C(H/3) and Wabs calculated iterativcly. See §3.2 for more details. 
In A. 

In units of lO'^K, calculated using empirical calibrations of lPilvuginI I2OOII) 

^ In units of lO'^K, Te[OII] calculated from empirical calibration and the relation given bv lGarnetj 1199211 
^ Radial velocity with respect to the center of Mkn 1087, in km s~^ 
^ In units of 10~^^ erg s~^ cm~^ 

" Estimated ages using our spectroscopical and photometric data, in Myr. See §4.1 for more details. 



smaller than the photometric uncertainties. The choice be- 
tween the two reddening curves implies only slight changes 
in the results and does not change our conclusions. The 
dat a in U. B . i?, J, H and Kg filters were corrected using 
the lRieke fc Lebofskvl |l98^ extinction correction law for 
Ajj, Ab, Ah, Aj, Ah and Aks- For the main body of Mkn 
1087 we have assumed an average extinction using the ob- 
jects C, #7, #8, #13, #14 and #15, that are the knots 
inside Mkn 1087 for which we have spectra. ISchleeel et all 
l|l998l) presented an all-sky map of infrared dust emission, 
determining a color excess of E{B — V) = 0.07 (assuming 
Rv = 3.1) at the Galactic longitude and latitude of Mkn 
1087. This value should be considered as a lower limit to 
the real one. We have assumed only this Galactic contri- 
bution to the reddening in the direction of Mkn 1087 for 
the more external objects (#1, #2, #3 and #12) and the 
galaxy KPG 103a because we have no determination of 
C(H/3) for it. The final E{B - V) adopted for each indi- 
vidual object is also shown in Tabled We have determined 
the errors for the photometry considering the FWHM of 
the PSF, sky level and flux calibration for each frame. The 
error in V — J color due to the differences in the shape 



and size between the optical and the NIR integrated areas 
is small when compared with the other uncertainties. 

3.2. Intermediate resolution spectra 

We have observed five slit positions in Mkn 1087 in order 
to obtain the spectra of the most interesting zones. In 
Figure 131 we show the four different position angles (P. A.) 
obtained using the INT over the B contour map, as well 
as the center of each extracted 1-D spectrum. All four 
slit positions cover the central zone, C, of Mkn 1087 and 
were combined to obtain an individual high signal-to-noise 
spectrum of this zone. The slit with P.A. 6° covers knots 
#1, #8 and the new bright emission area #13. The slit 
with P.A. 62° passes through knots #12 (that does not 
show nebular emission) and the new object #14. Knot 
#3 was observed with P.A. 120°, but it only reveals a 
faint nebular emission. The new companion galaxy and 
the intense knot #7 were observed using the slit position 
with P.A. 357°, that also shows the new nebular zone #15. 
The single slit position observed using WHT has a P.A. 
358° and is not represented in Figure El (it is similar to 
the P.A. 357° shown in the figure) and covers the center 
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Fig. 4. ISIS WHT spectra of the center of Mltn 1087, the intense knot #7 and its north companion object. The most important 
emission lines are labeled in the upper spectrum. All the spectra have been scaled down in flux in order to distinguish the faint 
lines. Note the high continuum level and the H I absorptions due to the underlying stellar population in the spectrum of the 
center of Mkn 1087. All spectra are redshift-corrected. 



of Mkn 1087, the companion galaxy and knot #7. The 
distance to the center of Mkn 1087 and the size of each 
extracted 1-D spectrum is indicated in Table |31 Note that 
we have observed six different areas at different radii along 
the main body of Mkn 1087, two of them (C and #7) 
observed both with INT and WHT. 

The wavelength and flux calibrated spectra of the cen- 
ter of Mkn 1087, the north companion and knot #7 are 
shown in Figure ^ The spectrum of the center of Mkn 
1087 shows a high continuum level and an intense under- 
lying stellar absorption in Hi Balmer lines. Absorption 
wings are clearly resolved for H/3, H7 and US. The high 
continuum level and these absorption features can be in- 
terpreted as the product of a substantial population of 
stars with ages larger than 10 Myr (see §4.1). The spec- 
tra of the other regions or knots show an almost absent 
continuum level with no evident absorption features, char- 
acteristic of a H II region-dominated zone. 

Fluxes and equivalent widths of the emission lines for 
each spectrum were measured using the SPLOT routine of 
the IRAF package. This task integrates all the flux in the 
line between two given limits and over a fitte d local contin- 
uum. We also used Starhnk DIPSO (,Howarth fc Miirravl 



software for the Gaussian fitting of profiles affected 
by line-blending (mainly, the [On] AA3726,3729 doublet in 
WHT spectra) and to obtain the best values for the fluxes 
of H I Balmer lines. For each single or multiple Gaussian 
fit, DIPSO gives the fit parameter (radial velocity cen- 
troid, Gaussian sigma, FWHM, etc) and their associated 
statistical errors. To correct the observed line intensities 
for interstellar reddening we have used the reddening con- 
stant, C(H/?), derived from the fluxes of Hi Balmer lines. 
However, we must also take into account the effect of the 
underlying stellar absorption in the hydrogen lines. Their 
equivalent widths, Wabs, are usually considered, but their 
values ar e uncertain beca use they also depend on the age 
of burst dOloffssonlhflfli . We have performed an itera- 
tive procedure to derive both, C(H/3) and Wabs, for each 
spectrum, assuming that the equivalent width of the ab- 
sorption lines is the same for all the Balmer lines and us- 
ing the relation given by iMazzarrella fc Boro nson' ('l993") 
to the absorption correction together with the .Whitforch 
(Eiii law: 



C{HP) 



1 



/(A) 



log 



WhS3 J 



TWW) 



(1 + ^) 



(1) 
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Table 4. Dereddened line intensity ratios with respect to I(H/?)=100 for the three brightest star-forming bursts 
observed in Mkn 1087: the center of Mkn 1087 (C), the north companion galaxy (N) and the intense burst #7. All 
have been observed both with 4.2m WHT and 2.5m INT. Correction for underlying stellar absorption in Hi Balmer 
lines has been also applied. 



Line 




f(A) 


C WHT 


C INT 


N WHT 


N INT 


#7 WHT 


#7 INT 


3726 


[0 III 


0.27 


206±1T' 


226±13'' 


155±12 


325±26° 


85.5±7.1 


212±20" 


3729 


[O III 

J 


0.27 






231±14 




106±8 




3869 


[Ne III] 


0.23 


8.3±2.9 


8.7±3.3 


43.4±6.9 


37.8±9.1 


9.9±3.6 


17.1±6.2 


3889 


He I-I-H8 


0.22 


13.8±4.3 


7.6±2.6 






15.0±4.3 


15.6±4.6 


3968 


[Ne III] 


0.21 


9.7±3.5 


3.54: 






11.0±3.6 


13.2±3.7 


4068 


[S II] 


0.19 


2.23: 












4076 


[S II] 

L J 


0.19 


1.13: 












4101 


m 


0.18 


25.9±4.3 


25.9±4.6 


25.9±4.8 


25.9±7.7 


25.9±4.4 


25.9±4.7 


4340 


H7 


0.135 


46.9±4.7 


46.9±5.2 


46.9±7.2 


46.9±8.8 


46.9±8.1 


46.9±9.8 


4363 


[0 ml 


0.13 


1.76: 










1.85: 


4471 


He I 


0.10 


2.14: 


3.90±1.5 






3.23: 




4861 


H/3 


0.00 


100±5 


100±6 


100±5 


100±7 


100±5 


100±6 


4959 


[o ml 


-0.02 


32.2±2.3 


33.9±2.7 


78.2±5.4 


124±7 


69.0±5.5 


57.7±6.1 


5007 


[o ml 

L J 


-0.03 


97.4±5.9 


109±7 


218±10 


375±19 


201±12 


161±10 


5016 


He I 


-0.03 


1.18: 












5200 


[NI] 


-0.05 




2.97±0.79 










5518 


[CI HI] 


-0.17 




0.94: 










5876 


He I 


-0.23 


10.3±1.2 


11.1±1.9 


12.6±2.2 


14.5±2.4 


10.3±1.3 


12.6±1.9 


6300 


[O I] 


-0.30 


6.73±0.60 


7.56±0.76 


19.1±3.0 


9.7±2.2 


4.55±0.79 


4.9±1.9 


6364 


[O I] 


-0.31 


1.83: 


1.77: 






1.41: 




6548 


[N II] 


-0.34 


31.4±2.7 


31.5±3.3 


8.0±1.7 


9.4±3.3 


14.2±1.9 


14.9±2.2 


6563 


Hq 


-0.34 


286±11 


286±12 


286±12 


286±14 


286±12 


286±14 


6584 


[N H] 


-0.34 


102.2±6.4 


98.1±8.8 


22.6±2.7 


22.8±4.1 


43.1±3.0 


46.5±5.4 


6678 


He I 


-0.35 


1.84: 






10.3±2.9 


2.09: 




6716 


[S II] 


-0.36 


42.9±3.8 


43.8±3.9 


48.7±4.9 


34.8±4.3 


31.0±3.4 


32.0±5.7 


6731 


[S II] 


-0.36 


36.0±2.9 


35.1±4.2 


37.4±4.1 


27.4±3.6 


23.7±2.2 


26.4±2.9 



" [O II] AA3726,3729 emission flux. 



for each detected hydrogen Balmer line. In this equa- 
tion, F{X) and /(A) are the observed and the theoretical 
fluxes (not affected by reddening or absorption), whereas 
Wabs , Wx and Whp are the equivalent widths of the un- 
derlying stellar absorption, the considered Balmer line and 
H/3, respectively. /(A) is the reddening curve normalized 
to H/3. We have considered the theoretical ratios between 
each pair of H l B almer lines expected for case B recom- 
bination given bv lBrocklehuriil Jl97lh and appropriated 
electron temperatures between 6000 and 10000 K and elec- 
tron densities around 100 cm~^. The C(H/3) and Wabs that 
provide the best match between the corrected and the the- 
oretical line ratios for each knot are shown in Table |21 We 
have assumed Wabs = 2 A in knots #13 and #15 because 
we only detected the two brightest Hi Balmer lines. For 
#1 and #3 we have not corrected for Wabs due to their 
faintness absorption. 

Finally, all the emis sion lines detect ed were corrected 
for reddening using the'Whi tfordI l)l958l) law and the cor- 
responding value of C(H/3). The H/3 line flux, F{H(3) (cor- 
rected by reddening and underlying stellar absorption) 
and the equivalent width of several lines (Ha, H/3, H7 and 
[Oiii]) are given in Table |3 while the reddening-corrected 



line intensity ratios relative to H/3 of each observed ob- 
ject are given in Tables^ and We have estimated the 
uncertai nties in the l i ne int ensities following the equation 
given bv lCastellanosI l)2000() . Colons indicate uncertainties 
of the order or greater than 40%. 

3.2.1. Physical conditions of the ionized gas 

The 1-D spectra were used to study the physical con- 
ditions and chemical abundances of the ionized gas. 
Unfortunately, the auroral [O ill] A4363 line is not detected 
or well measured, so we can not derive the electron tem- 
perature of the ionized gas. Although [Oiii] A4363 seems 
to be detected in the spectrum of the center of Mkn 1087, 
the broad absorption wing in H7 does not permit a proper 
measurement of this weak emission line. Consequently, the 
oxygen abundance for the different knots in Mkn 1087 was 
derived using the so-called empirical calibrations (see next 
section). The electron temperatures have been estimated 
from the T([Oiii]) and T([Oii]) pairs that repr oduce the 
total oxygen abundance obtained applying the IPilvueinl 
l|200ll) empirical method (see next section). As Pilyugin, 
we have assumed a two-zone scheme, the usual relation 
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Table 5. Dereddened line intensity ratios with respect to I(H/3)=100 for the faint knots observed in Mkn 1087. 
Correction for underlying stellar absorption in Hi Balmer lines was also applied. All spectra were taken using INT. 



Line 




#1 


#3 


#8 


#13 


#14 


#15 


O r-rc\ i—r \ TTl 

3727 [O IIJ 


0.27 


265: 


240: 


290±35 


252±75 


271±26 


286.39: 


3889 He I+H 8 


0.22 






19.5±6.4 








4340 H7 


0.135 






46.9±7.9 




47±14 




4861 H/3 


0.00 


100: 


100: 


100±9 


100±23 


100±12 


100: 


4959 [0 III] 


-0.02 






56.7±8.5 


55±14 


50.7±6.6 


34.7: 


5007 [0 III] 


-0.03 


136: 


99: 


165±14 


144±22 


173±14 


86.5: 


5876 Hel 


-0.23 










11.4±4.1 




6548 [N II] 


-0.34 






16.8±5.9 


26±10 


13.0±4.5 




6563 Hq 


-0.34 


286±46 


286±54 


286±23 


286±31 


286±10 


286±34 


6584 [N II] 


-0.34 


54±18 


51.5: 


60.0±9.1 


51.0±7.1 


59.7±7.2 


71±25 


6716 [S II] 


-0.36 






45.7±7.4 


46.7±7.6 


43.2±6.5 


22: 


6731 [S II] 


-0.36 






36.0±7.9 


36.0±5.8 


27.4±5.4 


57: 



Table 6. Chemical abundances of the knots. 



Knot 


C WHT 


C INT 


N WHT 


N INT 


#7 WHT 


#7 INT 


#1 


#3 


#8 


#13 


#14 


#15 


12+log 0/H° 


8.57 


8.53 


8.23 


8.24 


8.55 


8.53 


8.50: 


8.52: 


8.41 


8.47 


8.43 


8.44 


12+log O/H'' 


8.79 


8.78 


8.31 


8.32 


8.52 


8.54 


8.59 


8.58 


8.63 


8.57 


8.62 


8.68 


(0+++0+)/0 + 


1.55 


1.53 


1.53 


1.99 


2.19 


1.85 






1.57 


1.64 


1.61 


1.67 


12+log N/H 


7.76 


7.71 


6.77 


6.90 


7.48 


7.43 






7.33 


7.48 


7.32 


7.32 


- log N+/0+ 


0.81 


0.83 


1.46 


1.34 


1.08 


1.10 






1.08 


0.99 


1.11 


1.00 


12+log Ne++/H+ 


7.67 


7.60 


7.53 


7.42 


7.56 


7.79 














12+log S + /H+ 


6.54 


6.50 


6.24 


6.07 


6.31 


6.34 






6.38 


6.44 


6.33 


6.45 


12+log S + + /H+ 


6.42 
























12+log Hc+/H+'^ 


10.84 


10.88 


10.97 


11.03 


10.82 


10.94 










10.91 




[O/H] / [O/H] 0'' 


0.76 


0.69 


0.35 


0.35 


0.72 


0.69 


0.72: 


0.65: 


0.60 


0.71 


0.55 


0.56 



" Determined using the empirical calibration of IPilvugiiJ 120011) 
* Determined using the [N II] /Ha ratio IDenicolo et alj|200^ 

^ Using only the He I A5875 emission line^ 

Asuming 12+log[O/H]0 = 8.69+0.05 lAllende Prieto et alJl2001fl 



O/H = 0+/H+ + 0++/H+ and the linear relation be- 
tween T([Oiii]) an d T([O ii]) based on photoionization 
models obtained bv lCarnett (,1992,1 : 

r([0 ii]) = 0.7-T([O iii]) + 3000. (2) 

The empirical relation between t h e oxyg en abundance and 
the T([Oiii]) given bv IPilvueinl l)200l() (his equation 11) 
gives similar T[0 III] results than the ones obtained before. 
The finally adopted T[0 ill] and T[O i l] for each burst are 
shown in Table |31 FoUowing lPilvuginl l)200lj) . we could as- 
sume an uncertainty between 500 and 1000 K for T[0 ill]. 

We have used the [Six] AA 6717,6731 doublet to deter- 
mine the electron density, N^, of the ionized gas. Electron 
densities were calculated making use of the five-level pro- 
gram for the analysis of emission- Une nebulae that is in - 
cludcd in IRAF NEBULAR task ( Shaw fc Dufour"l 995|) . 
For zones jflA and #15, electron densities are below the 
low-density limit (<1GG cm^"^). We have no determination 
of A^e for #1 and #3. The final adopted values of Tg and 
Ne (or their upper limits) are compiled in Table 01 



3.2.2. Abundance Analysis 

We have derived the total oxygen abundance making use 
of empirical calibrations, base d on the r elatiy e intensi- 
ties of strong optical lines. D enicolo et al 1 (I2OO2I) used the 
[Nil] A6583 /Ha ratio in combination with photoioniza- 
tion models to obtain the calibration: 

12 + ,„g2.« + 0.73,ogEi4555^ (3) 

H Ha 

However, the most widely used empiri cal abundance cal i- 
brator is the i?23 parameter, defined as (|Pagel et al.ll980|) : 

_ /([O II] 3727) + /([O III] 4959 + 5007) 

= TW) ■ 

IPilvuginI ()2000() found that the previous calibrations us- 
ing this parameter had a systematic error depending on 
the hardness of the ionizing radiation, suggesting that the 
excitation parameter, P, is a good indicator of it. This 
parameter was defined as: 

/([O III] 4959-^5007) 1 . . 
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Fig. 5. a, b, c and d: Position-velocity diagrams for the four slit positions observed with INT, analysed in 1.6 arcsec bins. 
The horizontal bars represent the uncertainty of the Gaussian fitting for each point. The location and extension of the different 
galaxy members are also indicated. North is up in all the diagrams. 



performed a detailed analysis of the 
observational data combined with photoionization mod- 
els to obtain empirical calibrations for the oxygen abun- 
dance. Pilvugin 1,2000.,') proposes a linear fit that involves 
only the R-??. pa rameter (for 12+log(0/H)> 8.15), whereas 
jPilvuginl l)200lh is the following, more real and complex, 
calibration involving also the excitation parameter: 

(6) 



12 



, O 



i?23 + 54.2 + 59.45P - 



6.01 + 6.71P -I- 0.371P2 + 0.243i?23 ' 



(the so-called P-method), that can be used in moderately 
high-mctallicity H ii regions. 

In Table we show the oxygen abundances obtained 
for each observed knot using both empirical calibrations. 
We have assumed an error of ~ =^0-2 d e x for the oxygen 
abundances using iPenicolo et all l|2002l) . jPilvugini l|200lh 
estimates an uncertainty of ±0.10 dex for his calibration. 
We find that the O /H ratios obtained using the calibration 
of Dcnicolo ct al. ( 2002i) are systematically slight ly higher 
that t hose obtained using the calibrations of jPilvuginl 
ll200lh . Only in the brigh t burst #7 both c alibrations give 
a similar abundance. iDenicolo et all ((2002!) remarked that 



their calibration has a scatter that seems to depend on 
the degree of ionization. In consequence, we have adopted 
the oxygen abundances derived using the P-mcthod as the 
more appropiates ones. 

Note that the north companion object, N, has a clearly 
lower oxygen abundance [around 8.24, in units of 12-|-log 
(O/H)] compa red with the rest of the knots (between 
8.41 and 8.571. IVacca Contil ()lflfl2fl reported an oxygen 
abundance of 8.55 for the main body of Mkn 1087 assum- 
ing Te= 10000 K and 7Ve= 340 cm"^, value that coincides 
with the abundance we obtain for the central zone of the 
galaxy. The average value among all the observed knots 
inside Mkn 1087 is 12+log(0/H) = 8.49. 

In Table El we list the ionic and O abundances we have 
derived. All the ionic abundances have been calculated 
using the IRAF NEBULAR task from the intensity of 
collisionally excited lines, assuming the T([Oiii]) for the 
high ionization potential ions Ne++, and and 

T([Oil]) for the low ionization potential ions 0+, N+ and 
S+. 
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We have assu med the standard ionization correc- 
tion factor (icf) bv IPeimbert fc Coster o (1969), N/0 = 
N+ZO"*", to derive the nitrogen abundance of the objects. 
This correction factor is hsted in Table IHl Except for the 
central zone of Mkn 1087, only the S+/H+ ratio is derived, 
so we are not able to derive the total sulphur abundance. 
Furthermore, some contribution of S^"*" is also expected in 
the areas with a higher ionization degree. 

The Hc+/H+ ratio has been derived from the Hei 
A5875 line observed in each knot and using the predicted 
line emissivities calculated by .Smith et al,. (1996) . We 
have also correcte d for collisional contrib ution following 
the calculations bv lBeniamin et ;il ] l|2nn2ll . In Table El we 
present the Hc+/H+ ratios finally obtained. The results 
obtained from INT and WHT data for a given object are 
quite similar, despite the uncertainties. This indicates the 
good quality of the spectra and the reduction and analysis 
procedures. 

3.2.3. Kinematics of Mkn 1087 

The kinematics of the ionized gas were studied via the 
spatially resolved analysis of bright emission line profiles 
(Ha and II/3 for INT and WHT data, respectively) along 
each slit position. We have extracted zones of 4 pixels wide 
(1.6" long) covering all the extension of the line-emission 
zones in the four INT slit positions shown in Figure |21 
and zones of 6 pixels wide (1.2" long) for the WHT slit 
position. The Starlink DIPSO software was used to per- 
form a Gaussian fitting to the Ha and H/3 profiles. In 
figures and we show the position-velocity diagrams 
for the four INT slit positions and the single WHT slit 
position, respectively. All the velocities are referred to the 
mean heliocentric velocity of the center of Mkn 1087 (8336 
km s~^). The position of the different objects is also in- 
dicated on figures |31 and In Table |31 we show the mean 
radial velocities obtained for each zone, that span a nar- 
row interval of values with a maximum difference of about 
250 km s~i. 

All the position-velocity diagrams show a more or less 
clear rotation pattern in the inner part of Mkn 1087. Knot 
#7 participates of this rotation. However, the external 
zones have significant deviations from this behavior. Knots 
^1 and #3 are kinematically attached to the main body 
of the system despite their separation. We suggest that 
they are tidal dwarf galaxies (TDGs) formed from mate- 
rial stripped from Mkn 1087. Objects #13, #14 and #15 
seem to be kinematically attached to the rotation of Mkn 
1087. There is a small slightly sinusoidal pattern in the 
gas velocity around #13. This behavior, if real, could be 
produced by distortions associated w ith interaction effects 
i)Schweizeilll982HRubin et. al.lll990|) . Note that knot #1 
has a velocity similar to #13. The image of the system 
(see Figure^ reveals an optical bridge between both ob- 
jects, that furthermore have a very similar metallicity [12 
-|- log(0/H) ~ 8.50]. Figure El also shows some areas that 
can be interpreted as a tidal tail to knot #1. All these 
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Fig. 6. Position-velocity diagram for the single slit position 
observed in Mkn 1087 with WHT, analysed in 1.2 arcsec bins. 
The horizontal bars represent the uncertainty of the Gaussian 
fitting for each point. The location and extension of the differ- 
ent galaxy members are also indicated. North is up. 

evidences reinforce the idea that this structure (that has 
#1 and #13 at its extremes) could be actually a tidal tail. 

The HST ima ge of the center of Mkn 1087 
l|Malkan et alJll99^ shows a clear spiral pattern, imply- 
ing a dynamically mature system. It also reveals that its 
west side is nearer to us. From the position-velocity di- 
agrams of P.A. 357° and 358°, it is clear that the north 
side of the spiral pattern is the approaching one. If these 
two diagrams are interpreted as a circular rotation, then 
we can estimate the Keplerian mass of Mkn 1087. Taking 
into account the half of the maximum velocity difference 
(~125 km s"-'^) and the half of the spatial separation cor- 
responding to these maxima (~12" ~ 6500 pc), we obtain 
a mass of 2.37 x 10^° Mq for an inclination of 90°, as- 
suming circular orbits and Keplerian dynamics: 



M 



233 X d [pc] 



V [km s 



smi 



(7) 



A similar result is found using the position-velocity dia- 
gram of P.A. 6°. If we use P.A. 62° we estimate a veloc- 
ity difference of around 100 km s^^ in 8" (~ 4300 pc), 
obtaining a Keplerian mass of 1.0 x 10^*^ Mq, assum- 
ing an inclination of 90°. If we consider that the ellipti- 
cal shape of Mkn 1087 is only due to its inclination with 
respect the line of sight, we find that the proportion be- 
tween the mayor axis (that correspond with P.A. 357°) 
and the minor axis (between P.A. 62° and 120°) of the 
assumed disk is 5/3. This value corresponds with an incli- 
nation (defined as the angle betwe en the plane of the sky 
and th e plane of the galaxy) of 37° Gordon fc GottesmanI 
lll98ll) derived an inclination angle of 38°] . We now obtain 
a Keplerian mass of 6.53 x 10^" Mq for P.A. 357° (mayor 
axis) and 4.63 x 10^° Mq for P.A. 62° (minor axis). For 
the last one, we have also corrected the projected dis- 
tance, and it is the 70% of the mass obtained from the 
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position- velocity diagram of P. A. 357°. The average value 
between both estimations gives a Keplerian mass of Mkop 
= 5.6 X 10^° Mq for Mkn 1087. This value corresponds 
to a light to Keplerian mas s ratio of is/M^ep ~ 2.19. 
iGordon fc GottesmanI l)l98l[l . using the Arecibo antenna, 
estimated a total hydrogen mass of A/hi = 1.72 x 10^" 
M© and a total dynamical mass of A/dyn — 1-78 x 10"'^^ 
M© [fr om the velocity width of the H i profile and apply- 
ing the iBrandt | l|l96fl() rotation law] , that give a light to 
Hi mass ratio of Lb/Mhi ^1 .VI and a light to dynamical 
mass ratio of Ls/Mdyn ~ 0.69. Our Lb/Mj^cp ratio agrees 
with these estimates. Typical Lpi/ Myi^ values for spiral 
galaxies are between 11.2 and 8.9 ifBettoni et al.l 12003.) . 
Consequently, Mkn 1087 is a relatively gas-rich system. 



4. Discussion 

4.1. Ages of the bursts and stellar populations 



iKeell l)l988|) obtained position-velocity diagrams using 
various emission lines (Ha, H/3, [Nil], [Qui]), instruments 
and detectors for P.A. 57° (this slit position crosses knot 
#2) and for P.A. 68°. The last one is very similar to 
the position- velocity diagram we obtain for P . A. 62 °. The 
outer parts of the emission detected bv lKeell l)l988(l (that 
was chosen because it corresponds to one of the observed 
filaments) are redshifted by '^200 km with respect the 
inner part of Mkn 1087. He suggested that this feature in- 
dicates that this filament is infalling. However, as we will 
comment in §4.6, we think that the observed filaments in 
Mkn 1087 are tidal tails originated due to interaction with 
the nearby galaxy KPG 103a and the north companion 
object. 



The companion object has a radial velocity similar to 
the main body of Mkn 1087 (-1-117 km s"^), but it is kine- 
matically detached of the rotation pattern of the main 
galaxy. In Figure we can observe that the north com- 
panion object seems to have a sort of small rotation pat- 
tern, although its southern tip shows -if real- a deviation 
that could be interpreted as a feature of the interaction 
between it and Mkn 1087. Assuming Keplerian rotation 
and an inclination of 90°, we estimated a mass of 2.2 x 
10^ Mq for this object, just two orders of magnitude lower 
than the value estimated for Mkn 1087. Typical values of 
the mass for dwarf galaxies are between 10'' and 10^ Mq, 
so we can classify it as a dwarf object. These kinematic 
results, together with its low oxygen abundance ('^8.24) 
reinforce the evidence that it is an external galaxy and 
not a TDG. 



We remark the fact that knot #7 shows no important 
deviation of the main rotation pattern of Mkn 1087 in 
diagrams with P.A. 35 7° and 358°. This fact does not 
support the lKeell l)l988|) suggestion that #7 is consequence 
of the merger of a dwarf galaxy. In the best case, this 
presumed merger should be in a very advanced stage. Our 
data is more consistent with the interpretation of knot #7 
as a giant H ii region (or a complex of bright star-forming 
regions) located in one of the spiral arms of Mkn 1087. 



MEOO used lLeitherer fc Heckmanl l)l995|) models to derive 
the ages of the bursts using their optical photometric data. 
We have improved their results combining our reddening- 
corrected broad-band optical a nd near-infrare d photomet- 
ric data with STARBURST99 ( Leitherer et al.a99 91 mod- 
els (the improvement of Leitherer fc Heckman models). 
Spectral synthesis models with two different metallicities 
were chosen (Z/Zq = 1 and 0.4, the appropriate range of 
metallicities for the objects, see Table EJ, both assuming 
an instantaneous burst with a Salpeter IMF, a total mass 
of 10^ Mq and a 100 Mq upper stellar mass. The differ- 
ent color-color diagrams gave a similar interval of ages for 
all objects. In Figure [7^ the predictions of {V — R) ver- 
sus {U — B) are compared with our observational values. 
Knots with redder colors (#12, #3 and the bright galaxy 
KPG 103a) correspond to the older objects in the system 
(knot #11 should be included in this "red-group"). We 
have estimated an age of around 150 Myr for these three 
knots, values similar to those obtained by MEOO. 

A detailed analysis of the procedure used to derive 
the age is described in our study of th e Hickson Compact 



tne age is aescribea m our study oi tn e Miclison L^ompact 
Group 31 ijLopez-Sanchez et al. l2004j) . We have also used 
the W olf-R ayet starburst i nodels by ISchaerer fc Vaccal 
(Il998t) and IStasihska et alJ l)200l|) models to derive the 



age of the bursts using spectral features. A good agree- 
ment in all the age estimations has been obtained. We 
show the finally adopted value of the age of each spectro- 
scopically studied burst in TableO We have found that all 
the starbursts are between 6 and 9 Myr old, except knot 
#7 that seems to be younger (4.5 Myr). 

In FiguredD we plot our observational values of W(Ha) 
(obtained from the spectroscopic data) versus the color 
(U— B) (derived from the photometric data), and the pre- 
dictions of the STARBURST99 models. We find again a 
good agreement with the models except for knot #3, an 
older object that hosts a little star-formation activity. 

As we commented before, the three bright knots (Kl, 
K2 and K3) inside KPG 103a show blue colors [see their 
{U — B) colors in Table Although we have no spectra 
for them, they seem to be starburst zones in the nucleus 
of the ga laxy. W e have used the same STARBURST99 
(Leitherer et al.l ll999') models to estimate the age of the 
star-formation in them. As it can be seen in Figure 
the position of Kl, K2 and K3 are consistent with the 
models, and we estimate an age between 6 and 7 Myr for 
them. However, if the color excess, E{B — V), were higher 
than the one we have addopted (we only considered the 
Galactic contribution to the reddening for them, see §3.1), 
these knots will show bluer colors and, consequently, they 
will be a slightly younger. 

In Figure [S] we plot our observational values of 
ly (H/3) and the [O ill] A5007/H/3 emission hne flux showed 
in Tables El and 111 respectively, compared with the 
IStasihska et alJ l(200ll) models of H ii regions ionized by an 
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temporal evolution of the burst. The VF(Ha) of the north companion object is indicated with an horizontal arrow because there 
is not U ~ B color available for it. 



evolving starburst embedded in a gas cloud of the same 
metallicity. We indicate four models with two metalhci- 
ties, Z/Zq = 0.25 and Z/Zq = 1, changing the total mass 
and electronic density of the ionized gas. All the objects 
show a good fit with the predictions, except knots #13 
and #15. This fact can be due to a substantial population 
of old stars present in those objects, although it is diffi- 
cult to detect them in the spectra. However, as we saw 
before, their colors and M^(Ha) indicate that they host a 
starburst that are around 8 Myr old. 

The spectrum of the center of Mrk 1087 (see Figure^Ji) 
shows clear H i Balmer line absorptions, so we can estimate 
an effective age for the ste llar population in the nucleus 
using the STARBURST99 l)Leitherer et alJll999|) models. 
The best fit gives an age around 100 Myr, suggesting that 
star formation has been active in the nucleus of the galaxy 
ever since. 

In order to investigate the importance of older stellar 
populations, we have performed a simple analysis of the 
surface brightness profiles of the two main galaxies, as 
well as the north companion object and knot #3. We have 
taken concentric surfaces at different radii from the center 
of each system, and calculated the integrated flux inside 
each circle of area A using the relation: 



10 



fix = mx + 2.5 log ^, 



(8) 



to obtain the surface brightness, /i (in units of mag 
arcsec"^), of each circle A (in units of arcsec^) and mx 
is the magnitude in the filter X. In Figure O we show the 
surface brightness for B, V and R filters versus the size 
of the aperture (in arcsec) for Mkn 1087 and KPG 103a, 
whereas the surface brightness for the north companion 
and object #3 are shown in Figure ^| In these two fig- 
ures we also show the radial color profiles {B — R) and 
{B — V) (this last only in Figure EJ derived by a direct 
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Fig. 8. F([0 ml A 5007) versus W (H/3). Models by 
Istasiiiska et alJ i200lll . Tracks correspond to sequences of dif- 
ferent metallicities and electronic densities. Each symbol marks 
the position of the models at 1 Myr interval, starting in the 
upper-right corner of the diagram with an age of 1 Myr. 



subtraction of light profiles. The dotted horizontal line in- 
dicates the representative average color derived for each 
system (see Table 

We can observe that the surface brightness profile can 
be separated in two structures in all objects except #3, 
especially in the case of the B-band, indicating that they 
have a low surface brightness component underlying the 
starburst. In the case of the two brightest galaxies, we 
can attribute this underlying component to the galactic 
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Member 


U.C." 


Ms,o 


as (kpc) 




av (kpc) 


fJ-R,0 


an (kpc) 


Mkn 1087 


yes 


17.36 


1.63 


17.34 


1.72 


16.85 


1.55 


KPG 103a 


yes 


18.98 


1.64 


19.78 


2.44 


18.77 


1.76 


N. comp. 


yes 


20.91 


1.22 






20.89 


1.22 


#3 


no 


21.93 


1.42 






21.52 


1.41 



Underlying component. 




2 3 4 5 6 7 

R (arcsec) 
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Fig. 9. Surface brightness and B, V and R color profiles for 
Mkn 1087 (up) and KPG 103a (down). The line in the sur- 
face brightness diagrams is an exponential law fitting to the B 
profile. The dotted horizontal line in the color profile diagrams 
indicates the average color derived for each system. 

disk. We have performed an exponential law fitting to the 
profiles, following the expression: 



/ = /oexp(-Q!r), 



(9) 



that describes a typical disk structure: Iq is the central 
intensity and a is the scale length. The fitting structural 
parameters are indicated in Table [7| and the fit is plot- 
ted over each B profile with a straight line. It is not a 
surprise to find two structural components in Mkn 1087 
and KPG 103a: Figure [2| indicates the presence of a galac- 
tic disk with a possible spiral structure in both galax- 
ies. Furthermore, it is interesting to remark the noticeable 
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Fig. 10. Surface brightness and color profiles for the north 
companion object and knot #3. We follow the same notation 
than in Figure |^ Note that knot #3 has not an underlying 
component. 



variation in the {B — V) color along KPG 103a: it seems 
that the star-formation activity is very concentrated at 
its center, probably in the three inner knots (Kl, K2 and 
K3) analyzed in §3.1. Note that they show considerably 
bluer colors, (U — B) ~ —0.4, than the one derived for 
the complete galaxy, [U — B) = 0.17. We also find two 
structures in the north companion object (Figure in- 
dicating that an important old population is present and 
further supporting the idea that this object actually is an 
external galaxy nad not a tidal dwarf. However, we do 
not find two components in knot #:3. Is the contribution 
of the starburst dominating the underlying population? 
This problem of the "missing populations" in starbursts 
can b e solved using high spatial resolution photometry in 
NIR (|Noeske et al.ll2003|) . However, we do not think that 
it could be applied in #3: its low star-formation activity, 
colors and age suggest that most of the structure observed 
corresponds to an old stellar population. The weak contri- 
bution of the starburst is perhaps, in this particular case, 
diluted under the dominant old component. 

4.2. Star formation rates 
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and the classification of Mkn 1087 

Mkn 1087 and KPG 103a were detected jointly by IRAS 
(InfraRed Astronomical Satellite) at 25, 60 and 100 /im 
(source F04470+0314), with fluxes of /ss = 0.50 Jy, /go 
= 3.48 Jy and /loo = 4.77 Jy, respectively (IRAS Point 
Source Catalog, 1986). Using the Kennicutt (1998) cali- 
bration, we find that the star formation rate (SFR) de- 
rived for the system from infrared fluxes is SFRir=11.5 
Mq yr^^. It is a typical value for star formation galax- 
ies. We can also estimate the dust temperature. To, and 
the warm dus t ma ss, Mjj, using the expressions given by 
iBettoni et alJ l(200.11 and the /go and /loo fluxes. We ob- 
tain Td = 43.2 K and Md = 7.62 x 10^ Mq. Considering 
the luminosity of both galaxies, it corresponds to a mass 
to light ratio of log Md/ Lb = —4.36, a little lower than 
the expected values for spirals [between —3.74 and —4.07, 
iRettnni et ablEooi]. 

Radio emission from galaxies is aii extinction-free 
tracer of star formation. ICondon et al.l l)2002l) derived 
logLi.4GHz = 22.33 for Mkn 1087 using the NRAO VLA 
Sky Survey (NVSS). We can estimate the SFR from 
the 1.4 GHz luminosity using this data and the calibra- 
tion given bv lCondon et al" I liooi). We obtain a value of 
SFRi.4GHz(M > 5Mq) ~ 4.65 Mq yr^^, that could be 
considered as the rate of massive star formation during 
the past 10® yr. Assuming a Salpeter initial mass func- 
tion between 0.1 and 100 Mq, we can extrapolate it to 
include lower masses, finding SFRi.4GHz(Af > O.IAIq) ^ 
25.8 Mq yr"-'^, a relative high value compared with the 
SFR derived from FIR fluxes. This is a reasonable result 
because SFRi.4ghz measures the SFR of the most recent 
bursts of the star formation, while SFRir indicates the 
average SFR of a stellar population with larger average 
life times. 

The SFR can be also derived using the Ha emis- 
sion, applying the calibration given bv .Kennicutt U.998) . 
Although we do not have Ha imagery of the system, we 
can make an estimation of the Ha from our spectra con- 
sidering the slit size with respect to the total area of the 
bursts. We have estimated that the Ha flux (reddening 
corrected) we have measured in C is ~ 10 % of the total 
of the galaxy. This value corresponds to SFR^/q, ~ 4.5 Mq 
yr^^ for Mkn 1087. As we expect, this approximated value 
is lower than the ones derived from FIR and 1.4 GHz fluxes 
because it does not take into account the dust-obscured 
areas. The probable range for the real SFR in Mkn 1087 
is between 4.5 (from Ha flux) and 11.5 (from IR fluxes) 
Mq yr-1. 

Some authors, as lCondon et alJ l)2002() . classified Mkn 
1087 as Active Galactic Nucleus (AGN) but our data do 
not support it. The spectrum of the central zone of Mkn 
1087 (see Figure is similar to the ones of typical star- 
burst galaxies. The FWHM of the H/3 and [Oiii] A5007 
emission lines are 223 km s^^ and 228 km s^^ (corrected 
for instrumental broadening), respectively. The emission 
line widths in nuclear H ii regions and starburst galax- 
ies have values up to ~300 km s~^, whereas the typical 
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Fig. 11. Emission line ratios of knots in Mkn 1087. The solid 
lines give the limit for ionization by a zero age starburst, fol- 
lowing ,£o2ita^£t_aL, (.2000.) . 



FWHM me dian value range fo r AGNs is between 350 and 
550kms-i fv eron et al.'l997'). Furthermore, we can make 
use of the .Dopita et al.. (,2000) diagnostic diagrams based 
on emission line ratios to classify the excitation mechanism 
of the emission-line systems. In Figure these relations 
are plotted together with our observational data, being ev- 
ident that the knots studied here are consistent with the 
loci of typical H ii regions and not with AGNs. AGNs are 
clustered around log([0 ill] /H/3) ~ 1, log([Nii]/Ha) ^0 .1 

and log([Sii]/Ha) 0.2 llVeilleux fc Osterbrocklll987i) . 

However, Mkn 1087 could be classified as a Luminous 
Compact Blue Galaxy (LCBG). LCBGs are ~ L* {L* = 
1.0 X IO^'^Lq), blue, high surface brightness, high metal- 
licity, vigorous starbursting galaxies with an underlying 
older stellar population ( Guzman et al. 1998) . Specifically, 
a LCBG has {B — V) < 0.6, a surface brightness < 
21 mag arcsec"^ and Mb < -18.5. Mkn 1087 satisfies all 
these properties (for a radius of 5 arcsec, it has {B — V) 
= 0.11, a, fiB ~21 mag arcsec"^ and Mb ~ —22.3). 
LCBGs are quite common at intermediate redshifts, but 
by z their number density has decreased by a factor of 
ten and seem to have, on average, an order-of-magnitude 
smaller masses that typical galaxies today with similar lu- 
minosities (..Guzman et al. 2003) . The observed properties 
in LCB Gs at low z are si milar to the high z Lyman-break 
galaxies ( Erb et alJ2003(l . Consequently, LCBGs are espe- 
cially interesting for studies of galaxy evolution and for- 
mation because they could be the equivalent of the high z 
Lyman-break galaxies in the local universe and they have 
evolved more than other galaxy class. 



4.3. WR population 

We have not detected neither the blue Wolf-Rayet bump 
nor the nebular He ii A4686 emission line in the spectrum 
of the central zone of Mkn 1087. However , we have used 
the data obtained by IVaceli et al. 1 ^ll997^ to derive the 
WR population in the galaxy, assuming that all the flux 
came from the Heii A4686 emission line. These authors 
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estimated that F(H/3) = 7.81x10"^'^ erg s^^ cni"^ ^nd 
F(Heii A4686) = 8.36x10^1^ erg s^^ cm-^. if we consider 
that all the contribution of the Heii A4686 emission line 
comes from WNL we find around 7200 WNL stars assum- 
ing that -^(WNL 4686 ) = 1. 7x10^6 ergs"! for a WNL star 
JVacca fc Co nti 1992'). This v alue is somehow intriguing, 
because IVacca &: Contil l|l992l) reported an upper hmit of 
1600 to the total WNL presented in Mkn 1087, a factor 
4.5 smaller than our estir nate. Are rea ll y ther e Wolf-Rayet 
stars? The sht used by IVaceli et al.1 Jl997l) was 2. 4" x 
4.8", whereas the one used bv lVacca Sz Contil l)l992l) was 
1.5" wide, but they do not indicate its long. We used a 
1" X 5.4" aperture in our WHT observations. Perhaps, 
the WR sta r s are concentrated in a very small area and 
IVaceli et al ] (Il997l) were the only ones observing at the 
correct place. 

To derive the WR/(WR-|-0) ratio, we have to correct 
the contribution of the WR stars to the ionizing flux to 
obtain the total number of O stars, assuming that L(H/3) 
= 4.76x10^6 erg s'^ for a 07V star ifVacca fc Con ti 199 2) 
and r? = 07V/0 = 0.25 llSchaerer fc Vaccalll99^ . We de- 
rived around 47500 O stars in Mkn 1087, that implies 
a WR/(WR-t-0) ratio of 0.13. We have also used the 
calibration of the WR/ (WR +0) ratio using the flux of 
the blue WR bump given bv lSchaerer fc Vaccal l|l998|) in 
their study about evolutionary synthesis models for O and 
WR populations in young starburst (their equation 17). 
Assuming that all the blue WR bump flux comes from the 
Hen A4686 emission lin e and using the Hen A4686 flux 
from lVaceli et al.l l)l997|) . we estimate a WR/(WR-FO) ra- 
tio of 0.12, a value in excellent agreement with our pre- 
vious determination. However, these values represent an 
upper limit to the real one because IVaceli et al.1 ^^1997^ 
did not specify the stellar or nebular origin of the Hen 
A4686 emission line. 

The derived age (see §4.1 and TableOJ for the youngest 
bursts of star formation in Mkn 1087: knots C (the central 
zone) and #7, between 4 and 6 Myr, suggest that WR 
stars should be present in so young massive starbursts. 
Aperture effects seem to be playing an important role in 
the no-detection of WR features. 

4.4. The chemical composition of the objects. 

We have analyzed the radial dependence of the derived 
metallicity of the knots in Mkn 1087, which is shown 
in Figure ^] We have used the oxygen abundance ob- 
tained from the empirical calibration of Pilyugin (2QQ1J 
(see Table Ell, although this method gives a high uncer- 
tain oxygen abundance for knots #1 and #3 (we have 
assumed an error of ±0.20 dex for both, and ±0.10 dex 
for the rest of the objects, see §3.2.2.) 

We can distinguish two zones in Figure ^1 the first 
one is the galaxy disk itself, that extends to around 10 
kpc from the center of Mkn 1087. We find 5 different ob- 
jects in this area (the inner knots), two of them (C and 
#7) observed both with INT and WHT (we use the aver- 
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Fig. 12. Relation between metallicity and distance from the 
nucleus of Mkn 1087 for all the observed knots. We have con- 
sidered the oxygen abundances derived using the calibration 
bv .Pilvugin t.2001.1 . 

age value of both independent determinations). The object 
#15 lies just at the border of the disk. The second zone is 
the external area, where the north companion galaxy and 
knots #1 and #3 are located. 

We have explained above that knot #7 seems to be an 
intense star-formation zone off-center of Mkn 1087. It is 
resolved into nu merous individual kn ots in the HST im- 
age presented bv lMalkan et al.l l)l998l) . Figure IT^ suggests 
that the oxygen abundances of the center of the galaxy 
and #7 are very s i milar, despite the distance of 6 kpc be- 
tween them. iKeell l|l988! ) suggested that #7 is the conse- 
quence of a merger, but the position- velocity diagrams for 
P.A. 357° and 358° do not show any dynamical evidence 
of this. Furthermore, the individual knots that form ^1 in 
the HST image do not show a clear central concentration, 
supporting the idea that it is an interaction-triggered star- 
forming region and not the core of a stripped companion 
(| Walker et al.lil99 6). The oxygen abundance of the dif- 
ferent knots along the disk of Mkn 1087 is rather sim- 
ilar, but the central zone shows a slightly higher value 
[12+log(0/H) = 8.55, whereas the external zones show a 
values around 8.45]. On the other hand, the N+/0+ ratio 
is about 0.2—0.3 dex higher in the nucleus. These results 
suggest the possible presence of a rather weak abundance 
gradient along the disk of Mkn 1087 (about —0.013 dex 
kpc~^), although nothing definitive can be said taking into 
account the uncertainties. Typical values of oxygen gradi- 
ents in spiral galaxies are between —0.009 dex kpc'^ and 
—0.231 dex kpc~^ , with an average gradient of —0.06 dex 
kpc"i (jZaritski et al. 1994). 

The external knots #1 and #3 are located at 26 and 19 
kpc (respectively) from the center of Mkn 1087 and seem 
to show oxygen abundances similar to the ones obtained in 
the disk of the galaxy [12+log(0/H) - 8.50]. Taken at face 
value, this gives further evidence that knots #1 and #3 
are tidal dwarf galaxies. However, the high uncertainties in 
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their metallicities preclude any definitive conclusion from 
simple abundance considerations. 

From Figure IT^ it is clear that the north companion 
object has a significant lower 0/H ratio than the disk of 
Mkn 1087. The chemical difference is also evident in the 
N+/0+ ratio. is about —1.4 in the companion and 

between —0.8 and —1.1 in the rest of the knots. The value 
of —1.4 is of the ord er of that typical for dwarf galaxies 
ijGarnett et al.lll997|) . This is a further evidence that the 
north companion is an external pre-existing galaxy. 



4.5. The luminosity-metallicity relation 



iRicher fc McCalll lIlQQfj) obtained the following relation 
between absolute B magnitude and oxygen abundance for 
dwarf irregular galaxies (for A'Ib > ^18, the most common 
definition for a dwarf galaxy): 



12 + log ^ ^ 5.67 - 0.167 x Me 
H 



(10) 



that we have used for the objects in Mkn 1087. In 
Figure IK^I we show this relation (with an extrapolation 
to higher luminosities). Wc have assumed that Mkn 1087 
has the metallicity of its center, but obviously its position 
is under the [Richer fc McCalll l(l99^ relation because it is 
not a dwarf galaxy. 

In Figure El we can observe that the position of 
the north companion is consistent with the metallicity- 
luminosity relation for dwarf galaxies (in fact, it is the 
only object that is consistent with it). Therefore, the 
north companion should be interpreted as an independent 
nearby dwarf galaxy. The presence of the tidal dwarf #1 
and its associated bridge, which is almost aligned with the 
companion galaxy, suggests that the north companion is 
in interaction with Mkn 1087. 

iDuc fc Mirabell l)l998() found that dwarf galaxies with 
metallicities higher than the ones expected from the 
metallicity-luminosity relation (i.e, dwarf galaxies located 
at the upper-right zone in Figure IT^ are TDGs. Knots 
#1 and #3 ar e located in this z one o f the diagram and 
away from the iRicher fc McC^ l|l99!4l relation for dwarf 
irregular galaxies. This behavior is consistent with their 
possible TDG nature, in the sense that they have a much 
larger metallicity than it is expected for normal galaxies 
of the same luminosity. This evidence reinforces our hy- 
pothesis that both objects are tidal dwarf galaxies formed 
from material stripped by Mkn 1087, perhaps due to the 
interaction with the dwarf north companion object (in the 
case of #1 seems likely) and maybe also with the neigh- 
bour galaxy KPG 103a. 

Finally, the position of knots #7 in the metall icity- 
lumin osity relation does not support the claim by iKeell 
l)l988|) that this object is a merged external galaxy. If this 
is true, the position of the object should be more consis- 
tent with the relation. 
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Fig. 13. Luminosity-metallicity relation for the observed de- 
tached bursts in Mkn 1087 and knot #7. The solid line is th e 
relation for dwarf irregulars found bv lRlcher fc McCalJ (|l99i4 ). 
while the dashed line is an extrapolation of it. We have consid- 
ered the oxyge n abundances derived using the calibration by 
|Pilvuginl(|200ll) . 



4.6. Mkn 1087: A probable group in interaction. 

Sometimes, interaction features are found in group of 
galaxies, especially in Hickson Compact G roups [see, for 
example, llglesias-Paramo fc VilcheJ ll200lh: HCG 92: th e 
Stevhan's Quintet ("Yuu ct al."l997^ 'Sulentic et al.ll200ll) : 
HCG 54 Iverdes-Montcnegro et al. 2002): N GC 6845 
galaxy group llScott. Koribalski fc Jonesll2003(l : HCG 31 
llLonez- Sanchez et al.l2004^ ]. We consider that Mkn 1087, 
KPG 103a and the north companion object are, in fact, a 
group of galaxies. We might also include the tidal object 
#3 in this group. 

W e ha ve studied the selectio n criteria given bv lHicksonI 
l(l982l) and lHickson et sl\ to define a compact group 

of galaxies in the particular case of Mkn 1087. Our group 
satisfied the criteria of the number of members (at least, 
three galaxies that have similar radial velocities), isola- 
tion (there is not other galaxies with Amag <3 in a radii 
lower than 3 times the assumed size of the group) and 
compactness criterion (the mean surface brightness of the 
group, fiG, defined as the distribution of the fiux of all the 
members over the size of the group, is around 23.5 and, 
therefore, lower than 26). However, it fails in the relative 
magnitude criterion because the north companion galaxy 
differs in more than 3 magnitudes from the absolute mag- 
nitude of the brightest galaxies. Consequently, Mkn 1087 
and its surroundings can not be strictly classified as a 
compact group of galaxies attending to Hickson's criteria, 
which also are defined by the limitations of Hickson's sur- 
vey. What it is clear to us is that they actually form a 
group of galaxies with, at least, three members (four, if 
we also consider object #3). 

The knowledge of the dynamics and galaxy orienta- 
tions in the group offers some clues to understand what 
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kinds of interactions can lead to the observed structures. 
Mkn 1087 is a spiral galaxy with an inclination of ^ 37° 
with respect to the visual. Its west side is nearer to us and 
the north side is the approaching one (see §3.2.3). All the 
tails seem to be nearly radial to the spiral disk, but the 
streamers pointing toward objects #11 and #12 might be 
roughly polar, while the tail to knot #1 could be close to 
the disk plane. It is intriguing that intense burst #7 is 
also aligned in this direction. 

Knot #15, located at the beginning of the long tail to 
object #1, shows blue colors, [U — B) ^ —0.6, indicating 
star formation activity. The images (see Figures H and 
reveal small knots along this tail, that decrease in number 
in the direction to #1. They seem to be tiny star- forming 
regions too, because very weak Ha emission is detected in 
this zone of the WHT spectrum. Knot #1 appears to be an 
irregular object composed by diffuse little knots spanning 
about 4". We do not think it is a disrupted pre-existing 
object instead of a tidally stripped remnant because its 
metallicity is very similar to the spiral disk of Mkn 1087. 
Furthermore, if knot #1 were a disrupted low-mass ob- 
ject, its gravitational potential would be strong enough 
to modify consider ably the spiral pattern of Mkn 1087 
llWalker et al In "minor" mergers, the main galaxy 

develops a strong two-armed spiral pattern, whereas the 
nucleus of the low mass external satellite will finally be 
a small bulgelike entity. Ex amples of minor mergers are 
Arp 31 l|van Moorsel"l98?^. NGC 72 52 [the Atoms -for- 
Peace galaxy (Sch wcjzer 1982 : Hib bard et al.lll994|) ] or 
NGC 1097 fais-don fc Wallinll2003^ . All these systems 
show that the inner part of the main galaxy is very dis- 
turbed with ripples, loops and tails. None of these features 
are observed in our case: the HST image of Mkn 1087 re- 
veals a clearly and undisturbed spiral pattern, something 
that is also observed in the position- velocity diagrams (see 
Figures El and El) , and knot #1 is no t a bulgelike object. 
Moreover. »Mihos fc HernauistI (Jl994') proposed that these 
kind of mergers can lead to a significant inflow of disk gas 
and fuel a strong central starburst in the main galaxy, be- 
ing predominantly concentrated in the inner kiloparsec of 
the disk rather than throughout the disk or bulge. This 
is not observed in Mknl087; the star-formation activity is 
found all over the disk. 

The new data presented through our study give some 
clues to understand the evolution of the group. The num- 
ber and intensity of tails and bridges that are observed 
indicate that interaction processes have been important 
in the recent past of the system. The gravitational distur- 
bances seem to be the triggering mechanism of the star- 
formation bursts that nowadays are observed in the knots 
of Mkn 1087. We propose the following scenario to ex- 
plain the morphological features and chemical abundances 
present in the group: 

1. Three independent galaxies form the group: Mkn 1087, 
KPG 103a and the north companion galaxy, this last 
one is less massive than the others. Mkn 1087 is a spiral 



galaxy nearly side-on. KPG 103a is probably another 
spiral galaxy. 

2. It is possible that Mkn 1087 and KPG 103a suffered a 
nearby encounter that formed the tidal features ob- 
served towards knots #11 and #12 from material 
stripped from Mkn 1087. We think that it also origi- 
nated the tidal tail towards knot #3, that then formed 
a tidal dwarf galaxy with an important old population 
(although some star-formation activity also occurred 
because it shows weak emission features). KPG 103a 
also developed some type of tidal features that now are 
difficult to detect; they are barely observed at the east 
of KPG 103a, although it also shows a diffuse plume 
towards the west. There are blue knots (possibly, star 
forming bursts) in the center of the galaxy, suggesting 
that perhaps KPG 103a is still suffering the effect of 
tidal perturbations. 

3. Some time later, the north companion object proba- 
bly interacted with the north-east area of Mkn 1087. 
This encounter formed the tidal tails towards knot #2 
and the intense bridge of material towards knot #1, 
perhaps triggering the starburst in both galaxies. This 
close encounter had to be intense and relatively recent, 
because the structures that originated them are clearly 
detected, and were made from material stripped from 
the main body of Mkn 1087. Some star formation ac- 
tivity is also found in the tail to knot #1. The very 
weak plume observed at the south of the north com- 
panion galaxy, that points towards knot #1, reinforces 
this interacting scenario. 

Both interacting processes were probably not timely 
separated, i.e, perhaps the interaction between Mkn 
1087 and the north galaxy began before the interac- 
tion with KPG 103a end e d. -/V-Body sini ulations [i.e. 
iMihos fc H ernauisi^ (^1996^; IP^ibinskil ^2^] show that 
the elliptical galaxies found in most clusters and groups of 
galaxies have formed through the merger of several galax- 
ies, although galactic cannibalism of smaller galaxies only 
accounts for a small fraction of the accreted mass. Hence, 
the probable future evolution of the system will be the 
merging of all galaxies and, maybe, an elliptical galaxy will 
be formed. It would be desirable to perform detailed stud- 
ies of further LCBG in the Local Universe (as Mkn 1087). 
This could be essential to understanding the origin and 
evolution of this kind of galaxies in the earlier Universe. 
Perhaps, galaxy merging and interactions played a more 
important role at large z. 

Finally, it should be very interesting to obtain H i ra- 
dio map of the system to analyze the distribution of the 
neutral gas presents in Mkn 1087 and its surroundings (in 
a similar way to the studies b y iHibbard fc van GorkonJ 
(199^1 or lHibbard et all l(200lh . for example). This study 
will give a better understanding of the past and future 
evolution of the group galaxies, as well as it could allow 
to confirm definitely whether the tidal features observed in 
optical imagery are originated in the interactions between 
the different members of this group. 
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Table 8. B-magnitude (not corrected by reddening) and {B~ 
R) color of some of the background galaxies (bg) found in the 
area surrounding Mkn 1087. The objects belonging to the non- 
catalogued cluster (cl) show redder colors. 



Object 




{B~R) 


bgl 


18.49±0.08 


0.61±0.12 


bg2 


19.21±0.08 


0.36±0.12 


bg3 


19.43±0.08 


0.29±0.12 


bg4 


19.67±0.08 


0.70±0.12 


bg5 


20.6±0.1 


0.4±0.2 






0.4±0.3 


bg7 


21.6±0.2 


0.1±0.3 


bg8 


21.8±0.2 


0.6±0.3 


cll 


20.8±0.1 


1.1±0.2 


cl2 


21.8±0.2 


1.0±0.3 


cl3 


21.9±0.2 


1.8±0.3 


cl4 


22.0±0.2 


2.3±0.3 


cl5 


22.5±0.2 


2.7±0.3 


cl6 


23.2±0.3 


2.8±0.5 



4.7. The non-catalogued cluster of galaxies. 

The cluster of galaxies found at the upper right corner of 
the image in Figure^is not catalogued in the NASA/IPAC 
Extragalactic Database (NED). In Table|Hlwe show the B- 
magnitude and {B — R) color of some of the background 
galaxies (labeled with "bg") found in the deep i?-image 
and observed in the _B-image compared with the galax- 
ies belonging to the cluster (labeled with "cl"). We per- 
formed a simple analysis of the surface brightness profiles 
of some of these background galaxies to confirm that they 
are not star-like objects. Despite the uncertainties due to 
the faintness of the galaxies with ms > 20, the objects 
of this non-catalogued cluster show redder colors that the 
rest of the background galaxies. 

It is possible to determine the redshift to the cluster 
using photometric redshift codes, for ex ample, using the 
[Fernandez- Soto ct al. ( 199.8.) or HyperZ ijBolzonella et alJ 
|2000() codes, that find the redshift of each object using a 
fitting procedure that compares the observed magnitudes 
with the expected ones from standard spectral energy dis- 
tributions (SED). However, to apply these methods to our 
cluster, it should be necessary to get the magnitude in, at 
least, another broad-band filter. Nevertheless, using the 
code that is been developed by Juncosa fc Gutierrez [ some 
previous results can be found in l.Tuncosa et al.l 1 20041) ]. we 
estimate a redshift between 0.2 and 0.4 for it. The study of 
redshift clusters gives important clues about both galaxy 
and cluster formation and evolution. 



5. Conclusions 

We have used deep optical and near-infrared imagery and 
deep optical intermediate-resolution spectroscopy to an- 
alyze the morphology, colors, physical conditions, kine- 
matics and chemical abundances of the galaxy Mkn 1087 



and its surrounding objects. We have obtained the oxygen 
abundances for several star-forming bursts using empirical 
calibrations. We have obtained consistent age estimations 
for the knots observed in Mkn 1087 using different obser- 
vational indicators and theoretical models. 

We have derived a Keplerian mass of 5.6 x 10^" Mq 
for the main body of Mkn 1087. We have detected a new 
companion object at the north, which has an oxygen abun- 
dance and N/0 ratio lower that the rest of the knots, a 
Keplerian mass of ^ 2.2 x 10* Mq, an old underlying stel- 
lar population, and it is detached of the rotation pattern 
of the main galaxy. All these evidences suggest that the 
north companion object is an external dwarf galaxy. We 
suggest that some of the objects (#1, #2, #3, #11 and 
#12) are TDG candidates. 

Mkn 1087 is not powered by an AGN but it could be 
classified as a LCBG because of its color, absolute magni- 
tude and surface brightness. Wolf-Rayet features are not 
detected in the brightest bursts, perhaps due to aperture 
effects. Luminous knot #7 hosts a strong starburst that 
should not be the consequence of a merger, more likely 
it is a complex of star-forming regions located in a spiral 
arm of Mkn 1087. 

In conclusion, the complex geometry of the filamentary 
structure of Mkn 1087 and all the photometric, chemical 
and kinematical results can be explained assuming that it 
is in interaction with two external galaxies: 

1. the relatively bright KPG 103a, that could explain the 
bridges, the non-stellar objects located between both 
galaxies (#11 and #12), and the tidal dwarf galaxy 
#3; 

2. and the new dwarf north companion, that could orig- 
inate the tidal features towards knot #2, the bridge 
between Mkn 1087 and the knot #1, and produce the 
star formation triggering in the knots found in them. 

We consider that Mkn 1087 and its surroundings should 
be classified as a group of galaxies in interaction. 

Finally, we have also serendipitously detected a non- 
cataloged cluster of galaxies with central position around 
a = 04h 49min 35sec, S = +03° 21' 04" and probably at 
a redshift between 0.2 and 0.4. 
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